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Abstract 



Let J and J* be subsets of N such that 0,1 £ J and G J*. For infinitely many n, 
let k = (ki, . . . , k n ) be a vector of nonnegative integers whose sum M is even. We find 
an asymptotic expression for the number of multigraphs on the vertex set { 1, . . . , n} 
with degree sequence given by k, such that every loop has multiplicity in J* and every 
non-loop edge has multiplicity in J. Equivalently, these are symmetric integer matrices 
with values J* allowed on the diagonal and J off the diagonal. Our expression holds 
when the maximum degree fc max satisfies fc max = o(M 1 / 3 ). We prove this result using 
the switching method, building on an asymptotic enumeration of simple graphs with 
given degrees (McKay and Wormald, 1991). Our application of the switching method 
introduces a novel way of combining several different switching operations into a single 
computation. 

1 Introduction 

Let ki tn be a nonnegative integer for all pairs (i, n) of integers which satisfy 1 < i < n. Then 
for each n > 1, let k — k(n) = (&i, n , • • • , k n>n ). We usually write ki instead of k i>n . Define 
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M = Y^i=i hi- We assume that M is even for an infinite number of values of n, and tacitly 
restrict ourselves to such n. 

For subsets J, J* of N, define G(k, J, J*) to be the set of all multigraphs on the vertex set 
{1, . . . , n) with degree sequence given by k such that the multiplicity of every loop belongs 
to J* and the multiplicity of every non-loop edge belongs to J. Loops contribute 2 to the 
degree of their vertex. Also define Q{k) to be the set of multigraphs on {l,...,n} with 
degree sequence k (and no restrictions on the multiplicities). 

Equivalently we may think of Q(k,J,J*) as the set of all n x n symmetric matrices 
A = (dij) with all diagonal entries in J*, all off-diagonal entries in J and with 

2Qm + ajj = kj (1.1) 

for % — 1, 2, . . . , n. Note that diagonal entries are weighted by two in the row sum. 
Let /c max = max™ =1 hi. We find an asymptotic expression for 

G(k,j,r) = \g(k,j,r)\ 

that holds when k mSLX = o(M 1 / 3 ). For r = 1, 2, . . . let 



i=l 



where [a\ r = a(a — l)(a — 2) • • • (a — r + 1) denotes the falling factorial. Then M 1 = M and 
M r < A; max M r _i for all r > 2. 

For % > 0, define 

(l, ifieJ, fl, ifieJ*, 

Xi= < yi = < 

I 0, otherwise; I 0, otherwise. 

Our main result is the following. 

Theorem 1.1. Let J and J* be subsets of N such that 0, 1 G J and G J*. Suppose that 
n — > oo ; M — )> oo and A;^ ax = o(M). Then 



G{k, J, J*) 

Ml ( , n M 2 , 1N M| M 2 4 



(M/2)!2 M / 2 A; 1 !---A; ri ! 2 M v z 2y 2M 2 4M 5 

M 2 M 3 
~ 2M 4 



M 2 \ 
+ ( X3 _ X2 + (k 3 m jM)j . 
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For convenience, we restate the formula in the regular case. 
Corollary 1.2. Suppose that k = (k,k, . . . ,k) with kn even, such that k = o{n 1 / 2 ) as 



n — > oo. Then 



where 

Q(k,n) = \{k- l){(-l) X2 (k - 1) + 2(-l)") + ^(6x 2 - Qx 3 + 1). 

Setting yi = x 2 = x 3 = we recapture the asymptotic expression for sparse simple graphs 
with given degrees presented in [121 Theorem 5.2]. Similarly, setting y 1 = 1, x 2 = x 3 = we 
obtain the asymptotic enumeration by degree sequence of sparse graphs with loops allowed 
but no multiple edges: this is the second expression in [U Theorem 1.5]. 

We remark that the conditions G J* and 0, 1 G J in Theorem |1.1| can be replaced by 
weaker conditions: namely, that J* is non-empty and that J has at least two elements with 
the smallest two consecutive. Let s denote the least element of J* and let t, t + 1 be the 
smallest two elements of J. This case is reduced to ours by subtracting s from each diagonal 
element, t from each off-diagonal element, and 2s + (n — l)t from each ki, provided the new 



degree sequence satisfies the conditions of Theorem 1.1 



Theorem |1.1| is proved using the switching method, building on an asymptotic enumer- 
ation of simple graphs with given degrees |12j . Our application of the switching method 
introduces a novel way of combining several different switching operations into a single com- 
putation. 

The remainder of the paper is structured as follows. The history of this asymptotic 



enumeration problem is briefly reviewed in Section 1.1, then the new switching theorem is 
given in Section [2] In Section [3] we describe 15 types of switchings on multigraphs, which 
are used to show that certain multiplicities of edges or loops are rare. These switchings 



are analysed together in Subsection 3.14 using the new switching theorem. In Section [4] we 



complete the enumeration with the help of some calculations performed in 

Finally, in Section [5] we show that a naive argument leads to a formula for G(k, J, J* 



which differs asymptotically by a constant factor from the result of Theorem 1.1 in the 



regular case. The constant factor takes two different values, depending on whether 2 G J. 
1.1 History 

The earliest work on this problem was that of Read [14], who found exact and asymptotic 
formulae when k\ — ■ ■ ■ — k n — 3, for all four combinations of J* G {{0}, N} and J G 
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{{0, 1}, N}. The best result for J* = {0}, J = {0, 1} in the sparse range is that of McKay 
and Wormald [12] who treated /c max = o(M 1//3 ); see that paper for a survey of the many 
earlier results on the 0-1 case. In addition to J* = {0}, J = {0, 1}, Bender and Canfield [2] 
found the asymptotics for J* = {0}, J = N when fc max = 0(1). (Although that paper allows 
nonzero diagonal entries, they contribute singly to the row sums, not doubly as we have it.) 

In [TT], McKay and Wormald considered J* = {0}, J = {0,1} in the dense domain 
defined by min{k,n — k — 1} > cn/\ogn for c > | and \lq — k\ = 0(n l l 2+e ) for all i, where 
k is the average degree. McKay [8] found a better error term under the same conditions, 
while Barvinok and Hartigan allowed a wider range of degrees pQ . Greenhill and McKay jl] 
added the option of loops in both the sparse and dense ranges by considering J* = {0, 1}, 
J = {0, 1}. Finally, McLeod and McKay |9j analysed the case of J* = {0}, J = N when 
k\ = ■ ■ ■ = k n > cnj logn for c > \. 

For sparse rectangular matrices which are not necessarily symmetric, see Greenhill and 
McKay [3]. 



2 The switching theorem 

In order to bound the number of multigraphs having some unusual properties, we will apply 
the method of switchings. It will be necessary to apply several switching types, which we 
could analyse one at a time as in [3]. Instead we now prove a generalized switching theorem 
that allows us to analyse all of them at once. 

Theorem 2.1. Let T = r(V, E) be a directed multigraph, and let a : E — > M + be a positive 
weighting of the edges of T. Fix a non-empty finite set C , whose elements we will call colours, 
and let c : E — > C be an edge colouring of T (which need not be proper). For all v EV and 
c E C denote by r~(v) the set of edges of colour c entering v, and by -T c + (f ) the set of edges 
of colour c leaving v. 

We introduce the set of variables {N(v) : v £ V} U {s(e) : e £ E} and consider the 
following system of linear inequalities on these variables: 

N(v) > 0, (ye V) (2.1) 
a(e) > 0, (e £ E) (2.2) 
s(e) < N(v), (veV, ceC) (2.3) 

^(e)s(e) > N(v), (v £ V 7 c £ C, r+(v) ^ 0). (2.4) 

eer+(v) 
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Let C(v) = {c G C : r~(v) 7^ 0} be the set of colours entering v. For each v G V which is 
not a source, let X c (v),c G C(v) be positive numbers such that J2cec(v) ^c{v) < 1. For each 
vw G E, define a(vw) = a(vw) / \ c t VV] ){w) ■ Extend this function to any directed path P by 
defining a(P) = Yl ee pCn(e). 

Now suppose that X,Y, Z C V satisfy the following conditions: 

1. Z ^ 0, Z C X, and Y n 2" = 0; 

2. If v & V is a sink of T, or if a(vw) > 1 /or some t>w G £7, i/ien dG^. 

For any W, W C V, define V(W, W) to be the set of non-trivial directed paths in r 
which start in W , end in W , and have no internal vertices in Y U Z . Then every solution 



of (2.1) -(2.4) satisfies 



E N{v) < ^nrM(P) N(v), (2.5) 

where maximums over empty sets are taken to be 0. 

Proof. The case of one colour appears in [6, Theorem 3], apart from an inconsequential 
difference in the conditions on Z. We proceed by reducing the general problem to one in 
which there is only one colour, then modifying Z slightly so that [6, Theorem 3] applies. 

Define r~(v) and r + (v) to be the set of all edges (regardless of colour) entering v or 
leaving v, respectively. For each edge vw G E, define s(vw) = s(vw)\ c ( vw )(w). If we weight 



inequality (2.3) by X c (v) and sum over c G C(v), we obtain 



eer-(«) 



Similarly, if r + {v) 7^ 0, then we can sum (2.4) over those c G C with -T c + (f) 7^ to obtain 



a(e)s(e) > N(v), (v G V, r + (v) ^ 0). (2.7) 

eer+(u) 



Together with the nonnegativity of N(v) and s(e), we have equations of the form of (2.1) 



(2.4) with only one colour. 

Next remove from r all edges vw where v G Z. Doing so can only weaken the conditions, 



by decreasing the left hand side of some inequalities in (2.6) or removing some inequalities 



in (2.7). Note that none of the quantities in (2.5) are changed by removal of these edges 



Now Z satisfies the requirements of P, Theorem 3] , and applying that theorem to the system 



defined by (2.6), (2.7) and the nonnegativity of N(v) and s(e) gives (2.5) immediately. □ 
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We now describe how Theorem 2.1 can be used for counting. 

Suppose we have a finite set of "objects" partitioned into disjoint classes S(v), where 
v G V for some index set V. Define N(v) = \S(v)\ for each v G V. Also suppose that for 
each c6C we have a relation \P C between objects: to be precise, \P C is a multiset of ordered 
pairs (Q, R) of objects. (We call \P C a switching and usually define it by some operation that 
modifies Q to make R.) 

Now define an edge-coloured directed multigraph r = (V, E) with vertex set V, where 
r has a directed edge vw of colour c if and only if (Q,R) G ^ c for some Q G S(v) and 
R G S(w). (There is at most one edge of each colour between any pair of distinct vertices 
of r.) For each vw G E let s'(vw) = \{{Q,R) G \P C : Q G S(v),R G <S(w)}|, counting 
multiplicities, where c is the colour of vw. 

Fix v £ V and c G C such that -T 6 + (f) 7^ 0. Suppose that for any Q G 5(f) there are at 
least a c (v) > objects R with (Q,R) G !^ c , counting multiplicities. Then 

E s'(e)>a c (v)N(v). 
eer+(v) 

Similarly, for fixed v G V and c G C, suppose that for every R G S(w) there are at most 
b c {w) > objects Q with (Q, -R) G \P C , counting multiplicities. Then 

s '( e ) ^ b c (w)N(w). 

eer<r(w) 

Defining s(vw) = s'iyw) /b c ( vw ){w) anda(wwj) = b c ( vw ) (w) /a c ( vw ) (v) we obtain equations (2.1)- 



( 2.4[ ). Theorem 2.1 can thus be used to bound the relative values of J2 v ex\^( v )\ anc ^ 



YIv&y 1^(^)1 ^ there is some Z for which X,Y,Z satisfy the requirements of the lemma. 
A common choice (which we will use here) is X = V, as then we bound the fraction of all 
objects which lie in {J vGY S( v )- 



3 Switchings on multigraphs 

Define 

Ni = max{ [log M] , [480M 2 /M] }, 

N 2 = max{ [log M] , [240M 2 2 /M 2 ] } , (3.1) 

N 3 = max{ [log M] , [240M 3 2 /M 3 ] }. 

(We have not attempted to optimise constants.) Given a multigraph Q, let £d(Q) denote 
the number of loops with multiplicity D and let cd{Q) denote the number of non-loop edges 
with multiplicity D, for D > 1. 
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Let 

Go = g(k, J U {4, 5, 6, . . .}, J* U {2, 3, 4, . . .}) 

be the set of all multigraphs with degree sequence k, and allowing all multiplicities except for 
those in {1} — J* on loops and those in {2, 3} — J on non-loops. Note that Q(k, J, J*) C Go- 
We also define the subsets Y, Z of Go by 

Y = Go - {Q e Go I £d(Q) = for D > 2, e D {Q) = for D > 4, 

h(Q) < iVi, ea(Q) < iV 2 , e 3 (Q) < iV 3 }, (3.2) 
£ = {Q E Qo | £ D (Q) = for D > 2, e D (Q) = for D > 4, 

< [iVi/21, e 2 (Q) < [iV 2 /2l, e 3 (Q) < [iV 3 /2l}. (3.3) 

We define 15 coloured switchings which act to reduce the number of loops or edges with 
high multiplicities, moving in steps from Y towards Z. These switchings are defined below, 
together with a description of when each should be used. Indeed, any given switching will 
only be used on multigraphs in Gq for which none of the switchings with a lower-labelled 
colour are applicable. An important property of all the switchings is that they do not create 
simple loops (with multiplicity 1) and they do not create non-loop edges of multiplicity 
2 or 3, as these may not be allowed for multigraphs in Go- 

For each switching colour c and multigraphs Q, R we will define the following parameters: 

• a>c{Q) is a lower bound on the number of ways in which a switching of colour c can be 
applied to Q; 

• b c (R) is an upper bound on the number of ways in which a switching of colour c can 
produce R; 

• if R can be obtained from Q by performing a switching of colour c, then we let a(Q,R) 
be an upper bound on b c (R)/a c (Q). (Note that the colour c is determined by the 
edge QR.) 



Finally in Section 3.14 we will analyse all the switchings at once by applying Theorem 2.1 



3.1 Switchings of colour 1 

A switching of colour 1 is used to reduce the number of loops of multiplicity equal to 2 or 
multiplicity at least 4. It is are applied to multigraphs Q e Go with L(Q) > [3M 1 / 2 ] , where 

L(Q)=£ 2 (Q) + J2^d(Q) (3.4) 

D>4 
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is the number of loops in Q with multiplicity 2 or multiplicity at least 4. The switching is 
described by the sequence (t>i, v 2 ) of distinct vertices such that there is a loop of multiplicity 
D\ at v\ and a loop of multiplicity D 2 at v 2 , with Di, D 2 G {2} U {4, 5, . . .}. 

Let m be the multiplicity of the edge from v\ to v 2 in Q (which may equal zero). The 
switching reduces the multiplicity of both loops by 2, and increases the multiplicity of the 
edge {vi,v 2 } by 4. This operation is depicted in Figure [TJ 

D x D 2 D x -2 D 2 -2 

Figure 1: A switching of colour 1 



O m + 4 O 



Suppose that i? can be produced from Q using a switching of colour 1. The number of 
ways to perform a switching of colour 1 in Q is exactly 

[L(Q)h > 9M 

as there is no restriction on the value of m. The number of ways that R can be produced 
using a switching of colour 1 is bounded above by Mj 4, which is a bound on the number of 
ways to choose an oriented edge with multiplicity at least 4. Hence we can set ai(Q) = 9M 
and bi(R) = M/4, giving 

a(Q,R) = ± 

3.2 Switchings of colour 2 

A switching of colour 2 is used to reduce the number of loops of multiplicity three to below 
[3M 1 / 2 ] . It is performed for multigraphs Q for which switchings of colour 1 do not apply and 
such that 1?,{Q) > [3M 1 ' 2 ]. The switching is described by the sequence (^1,^2) of distinct 
vertices such that there is a loop of multiplicity 3 at v\ and at v 2 . 

Let m be the multiplicity of the edge from V\ to v 2 in Q (which may equal zero). The 
switching removes the loop at v\ and at v 2 , and increases the multiplicity of the edge {vi, v 2 } 
by 6, as illustrated in Figure [2} 

Suppose that the multigraph R can be produced from Q using a switching of colour 2. 
The number of ways to perform a switching of colour 2 in Q is exactly [-^(Q)^ > 9M, and 
the number of ways that R can be produced using a switching of colour 2 is at most M/6 
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Figure 2: A switching of colour 2 



(since an oriented edge of multiplicity at least 6 determines the reverse operation, from R 
to Q). Hence we can take a 2 (Q) = 9M and b 2 (R) — M/6, leading to 

a(Q,R) = ± 

3.3 Switchings of colour 3 

A switching of colour 3 is used to reduce the number of loops of multiplicity 1 to below 
|~M 1//2 ] . It is applied to multigraphs Q for which switchings of colour 1 and 2 do not apply 
and such that i\{Q) > \M . The switching is defined by the sequence (v±, v 2 , V3) of distinct 
vertices such that each of v%, v 2 , v 3 has a simple loop in Q (of multiplicity 1), and none of 
the edges {vi,v 2 }, {^1,^3}, {^2,^3} ar e present in Q. The switching removes these three 
loops and joins the three vertices pairwise by simple edges (of multiplicity 1), as illustrated 
in Figure |3j 



1 




Figure 3: A switching of colour 3 



Suppose that the multigraph R can be produced from Q using a switching of colour 3. 
The number of ways to perform a switching of colour 3 in Q is at least 

[*i(Q)] 3 -0(Wi(Q) 2 ) =4(Q) 3 (i-o(i)) > \M'i\ 

while the number of ways that R can be produced using a switching of colour 3 is at most 
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k max M. Therefore we can take a 3 (Q) = \M' i l 2 and 6 3 (R) = £; max M, leading to 

o(l). 



//n r>\ ^^max 

a(Q,R) 



M 1 ' 2 



3.4 Switchings of colour 4 



A switching of colour 4 is used to reduce the number of non-loop edges of multiplicity greater 

1 /2 

than max{4, |~/c m ax]}. It is applied to multigraphs Q for which switchings of colours 1,2,3 
do not apply and such that E + (Q) > [4fcma X M 1 / 2 ] , where 



E+(Q) = e ^Q)- 

D=max{4, \kU^}+l 



(3.5) 



A switching of colour 4 in Q is described by a sequence [y i, W\, t>2, W2) of distinct vertices 
such that 

• the multiplicity of edge (v\,w\) in Q is D\ and the multiplicity of (v2,w 2 ) in Q is D 2 , 

1 /2 

where D 1 ,D 2 > max{4, |/c ma xl}, 

• the edges {^1,^2} and {^1,^2} have multiplicity zero in Q (they are non-edges). 

The switching reduces the multiplicity of edges (v\,wi) and (^2,^2) by one and gives mul- 
tiplicity 1 to edges {^1,^2} and {wi,W2}, as shown in Figure |ij 



D 1 -l 



Do - 1 



Figure 4: A switching of colour 4 or 5 



Now suppose that the multigraph R can be produced from Q using a switching of colour 4. 
The number of switchings of colour 4 that can be performed in Q is bounded below by 

4[£ + (Q)] 2 - 0(k 2 max E + (Q)) = AE + (Q) 2 (1 - o(l)) > 64 A; max M(l - o(l)) > 60k m3X M, 
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and the number of switchings of colour 4 that can produce R is at most 



max{4, [/Cmax]} 2 

Hence we can take a±{Q) = 60 A; max M and b A (R) = 2/c max M, leading to 

a(Q,R) = ± 

3.5 Switchings of colour 5 

A switching of colour 5 is used to reduce the number of non-loop edges of multiplicity at 

1 /2 

least 5 and at most [fcmax] • For a multigraph Q E Go, let 

I7c 1/2 1 
I "'max I 

E-(Q)= *d{Q)- 

D=5 

Switchings of colour 5 are applied to multigraphs Q for which switchings of colours 1, . . . , 4 
do not apply and such that E~(Q) > \3k max M 1 ^ 2 ~\ . 

A switching of colour 5 is described by a sequence (v%, Wi, v 2 , w 2 ) of distinct vertices such 
that 

• the multiplicity of W\} in Q is D 1 and the multiplicity of {v 2 , w 2 } in Q is D 2 , where 

D 1 ,D 2 e{S,6,...,\k 1 JL~\}, 

• the multiplicity of {v\,v 2 } and {wi,w 2 } in Q is zero (these are non-edges). 

This switching is also illustrated by Figure [4j but with different conditions on Di,D 2 as 
above. 

Suppose that R can be produced from Q using a switching of colour 5. The number of 
ways that a switching of colour 5 can be performed in Q is bounded below by 

4[iT(Q)] a - 0{kl^E-{Q)) = AE-{QY(l - o(l)) > SOA^M, 

and the number of ways to produce R using a switching of colour 5 is at most 

Hence we can let a 5 (Q) = 30A^ ax M and b 5 (R) = k^^M, leading to 
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3.6 Switchings of colour 6, 7, 8 



A switching of colour 4 + j is used to reduce the number of edges of multiplicity j to below 
|~M 5//6 ], for j = 2,3,4. Switchings of colour 4 + j are applied to multigraphs Q for which 
switchings of colours 1, . . . , 3 + j do not apply and such that 

e j {Q)> [M 5 / 6 ]. 

Given such a multigraph Q, a switching of colour 4 + j is defined by a sequence 

(V 1 ,W 1 ,V 2 ,W 2 ,.. -,Vj,Wj) 

of distinct vertices such that 

• (v r , w r ) is an edge of multiplicity j in Q for r = 1, . . . , j, 

• every edge {v r ,w s } with 1 < r ^ s < j has multiplicity in Q (it is a non-edge). 

The switching deletes these j edges of multiplicity j, and inserts a complete bipartite graph 
Kjj on {vi, . . . ,Vj} U {w\, . . . ,Wj}, with all j 2 new edges simple (that is, multiplicity 1). 
This operation is illustrated in the Figure [5] for the case j = 4. 




Figure 5: A switching of colour 8 



Suppose that the multigraph R can be obtained from Q using a switching of colour 4 + j, 
for j e {2, 3, 4}. The number of ways in which a switching of colour 4 + j can be performed 
in Q is bounded below by 

while the number of switchings of colour 4 + j which produce R is at most kH~ 2 M. Hence 
we can define ci4_|_j 

(Q) = iM 5 ^ 6 and & 4 +i(#) = k%~ 2 M for j G {2,3,4}. Since j > 2, this 

leads to 
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Before describing more coloured switchings, we prove a useful fact. 

Lemma 3.1. Suppose that Q G Go is such that no switching of colour 1 to 8 applies to Q. 
Then ei(Q) = {\-o{l))M. 

Proof. Since no switching of colour 1, 2 or 3 applies we have 

J2d£d(Q) = 0(k^M 1 ' 2 ) = 0(M 5 / 6 ), 

D>1 

and since no switching of colour 4-8 applies we know that 

De D{Q) < 3[M 5 / 6 1 + k max E+(Q) + E-(Q) 

D>2 

= 0(M 5 / 6 + ^ a 2 x M 1 / 2 ) 
= o{M). 

The result follows. □ 
3.7 Switchings of colour 9 

Switchings of colour 9 reduce the number of loops of multiplicity 2 or multiplicity at least 4, 
until this number is zero. They are applied to multigraphs Q for which the switchings of 
colours 1, . . . , 8 do not apply and such that L(Q) > 1. 

Let Q be such a multigraph. A switching of colour 9 in Q is defined by a sequence 
(vq,Vi,wi,V2,W2) of distinct vertices such that 

• there is a loop at vertex vq in Q with multiplicity D, where D = 2 or D > 4, 

• {vi,wi} and {t^u^} are simple edges (of multiplicity 1), 

• {^0) w i}) {^0)^2}, {^o^i}) {^0,^2} are ah non-edges in Q (with multiplicity zero). 

The switching reduces the multiplicity of the loop at Vq to D — 2, removes the two simple 
edges {vi,Wi} and {^2,^2} and inserts the four simple edges {v ,Vi}, {^0,^2}, {v Q ,Wi}, 
{^0,^2}, as shown in Figure [6j 

Suppose that the multigraph R can be produced from Q by a switching of colour 9. The 
number of ways to perform a switching of colour 9 in Q is at least 

L(Q)([2 ei (Q)] 2 -CK*L*ei(Q))) > \ L(Q)M 2 , 
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using Lemma |3~T , while the number of ways to produce R using a switching of colour 9 is 
at most M4. (We ignore the presence of the loop at vq in this upper bound, since no such 
loop exists when D = 2.) 

Hence we can let ag(Q) = \L{Q)M 2 and bg(R) = M 4 , leading to 



max 



L(Q) M 2 \ M 
since L(Q) > 1. 



3.8 Switchings of colour 10 

Switchings of colour 10 reduce the number of loops of multiplicity 3, until this number is 
zero. They are applied to multigraphs Q such that switchings of colours 1, ... ,9 do not 
apply and such that CziQ) > 1- 

Let Q be such a multigraph. A switching of colour 10 is defined by a sequence of distinct 
vertices (vo,vi,wi,V2,W2,Va,W3) such that 

• there is a loop of multiplicity 3 at vq in Q, 

• there is a simple edge {vj,Wj} in Q, for j = 1, 2, 3, 

• the edges {vq, Vj} and {vo, Wj} all have multiplicity in Q, for j = 1, 2, 3. 

The switching removes the loop of multiplicity 3 and the simple edges {vj, Wj} for j = 1, 2, 3, 
and adds the six simple edges {{v ,Vj}, {v ,Wj} \ j = 1,2,3}, as shown in Figure [7j 

Now let R be a multigraph which can be formed from Q by a switching of colour 10. The 
number of ways that a switching of colour 10 can be performed in Q is at least 

£ 3 (Q) ([2 ei (Q)] 3 -0(A4 ax ei(Q) 2 )) > \HQ)M\ 
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Figure 7: A switching of colour 10 



using Lemma [37TJ and the number of ways that R can be produced using a switching of colour 

(Q) M 3 and b 10 (R) = M 6 . This 



10 is bounded above by M 6 . Hence we can take a w (Q) 
leads to 

2M 6 



2*3 



a(Q,#) 



A- 5 

q l "'max 



We do not tackle single loops yet. Instead, the next two switchings reduce non-loop edges 
of high multiplicity down to zero. 



3.9 Switchings of colour 11 

Switchings of colour 11 reduce the number of non-loop edges with multiplicity 4 or multiplic- 
ity at least 7, until this number is zero. Switchings of colour 11 are applied to multigraphs 
Q for which switchings of colours 1, . . . , 10 do not apply and E(Q) > 1, where 

= e 4 (Q) + $> D (Q). 

D>7 

Let Q be such a multigraph. A switching of colour 11 in Q is defined by a sequence 

(v , w , Vi,Wi, t> 2 , w 2 , v 3 , w 3 ) 

of distinct vertices such that 

• the edge {v , w } has multiplicity D in Q, where D = 4 or D > 7, 

• edges {vj,Wj} are simple edges in Q (with multiplicity 1) for j = 1,2,3, 

• the edges {v ,Vj} and {w ,Wj} are non-edges in Q for j = 1,2,3 (with multiplicity 
zero) . 
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The switching reduces the multiplicity of {vq,w } to D — 3, removes the edges {vj,Wj} for 
j = 1,2,3 and adds the simple edges {vo,Vj} and {wo,Wj} for j = 1,2,3. This operation is 
illustrated in Figure [HJ 




Suppose that the multigraph R can be obtained from Q by performing a switching of 
colour 11. The number of ways to perform a switching of colour 11 in Q is at least 

2E(Q) ([2ei(Q)] 3 - 0(k£ tM *e 1 (Q) a )) > E(Q) M 3 , 



using Lemma 3.1, and the number of ways to produce R using a switching of colour 11 is at 



most & max M 4 . Hence we can let a n (Q) = E(Q) M 3 and bn(R) = /c max M 4 . This leads to 



P Ma / k & \ f h 3 



t /~\ T3\ max * s~\ I max \ 



max 



E(Q)M 3 \ M 2 J \ M 
3.10 Switchings of colour 12 

Switchings of colour 12 reduce the number of non-loop edges with multiplicity 5 or 6, until 
this number is zero. They are applied to multigraphs Q such that switchings of colours 
1, . . . , 11 do not apply and e$(Q) + e§(Q) > 1. 

Let Q be such a multigraph. Then a switching of colour 12 is defined by a sequence of 
distinct vertices 

(v Q , w , hi, wt, v 2 , w 2 , v 3 , w 3 , v 4 , w 4 , v 5 , w 5 ) 

such that 

• the edge {vo,wo} has multiplicity D in Q, where D e {5,6}, 

• each edge {vj, Wj} is a simple edge in Q, with multiplicity 1, for j = 1, 2, 3, 4, 5, 
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• the edges {v , Vj} and {w , Wj} all have multiplicity in Q for j = 1, 2, 3, 4, 5 (that is, 
they are non-edges). 

The switching reduces the multiplicity of the edge {v ,w } to D — 5, removes the edges 
{v j, Wj}, j = 1, 2, 3, 4, 5, and inserts the simple edges {vo, Vj}, {vo, Wj} for j = 1, 2, 3, 4, 5, as 
shown in Figure |9j 

D D-5 




Figure 9: A switching of colour 12 



Suppose that the multigraph R can be formed from Q by a switching of colour 12. The 
number of ways that a switching of colour 12 can be performed in Q is at least 

2(e 5 (Q) + e 6 (Q)) ([2ei(Q)] 5 - 0(*4 ax ei(Q) 4 )) > (e 5 (Q) + e 6 (Q)) M 5 , 



using Lemma |3.1[ and the number of ways that R can be produced by a switching of colour 
12 is at most Mf . Hence we can set a 12 (Q) = (e 5 (Q) + e 6 (Q))M 5 and b 12 (R) = Mf . This 
gives 

3.11 Switchings of colour 13 

Switchings of colour 13 are used to reduce the number of simple loops until this number is 
at most [A^i/2] . They are applied to multigraphs Q for which switchings of colours 1, . . . , 12 
do not apply and £i(Q) > [A^i/2]. 
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Let Q be such a multigraph. Then a switching of colour 13 is defined by a sequence of 
distinct vertices (vq, vi, t> 2 ) such that there is a simple loop on t> in Q, the edge {v±, t> 2 } is a 
simple edge in Q, and the edges {vo,vi}, {^0,^2} are both absent in Q (that is, they have 
multiplicity zero). The switching removes the simple loop at vq and the simple edge {v±, t> 2 } 



and inserts the two simple edges {t> , t>i}, {t>o, f 2 }, as shown in Figure 10 




Figure 10: A switching of colour 13 



Suppose that the multigraph R can be obtained from Q by performing a switching of 
colour 13. The number of ways that a switching of colour 13 can be performed in Q is at 
least 

e l (Q){2e 1 (Q)-0(k 2 m J)>le 1 (Q)M, 



using Lemma |3.1[ and the number of ways that a switching of colour 13 can produce R is 
at most M 2 . Therefore we can define a>i 3 (Q) = Mx(Q)M and 613 (R) = M 2 , leading to 



2M 2 4M 2 1 
a[Ql ' ~ h(Q)M ~ iV\M " 120 : 



using the definition of N- 



3.12 Switchings of colour 14 

Switchings of colour 14 reduce the number of non-loop edges with multiplicity 2 until this 
number is at most [A^ 2 /2]. This switching is applied to multigraphs Q such that switchings 
of colours 1, . . . , 13 do not apply and e 2 (Q) > [iV 2 /2]. 

Let Q be such a multigraph. A switching of colour 14 in Q is described by a sequence 
(vq, w , vi, wi, v 2 , w 2 ) of distinct vertices such that 

• {^o, wo} is an edge of multiplicity 2 in Q, 

• {vi,wi} and {f 2 ,w 2 } are simple edges in Q (with multiplicity 1), 

• none of the edges {vo,Vj} or {wo,Wj} are present in Q, for j = 1,2 (these are all 
non-edges, with multiplicity zero). 
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The switching removes these 3 edges and replaces them with two copies of Ki 2} one on 
{vq,Vi,V2} centred at vq and the other on {wq,wi,W2} centred at wq. This operation is 



shown in Figure 11 



Figure 11: A switching of colour 14 



Suppose that the multigraph R can be produced from Q using a switching of colour 14. 
The number of ways to perform a switching of colour 14 in Q is at least 

2e a (Q)([2ei(Q)] 2 - 0(£; r 2 nax ei(Q))) > e 2 (Q)M 2 , 



using Lemma |3.1[ and the number of ways that R can be produced using a switching of 
colour 14 
leading to 



colour 14 is at most Mf. Therefore we can take a u (Q) = e 2 (Q)M 2 and bi^R) = Mf , 



Mf 2M| 1 

tt(g ' } ~ e 2 (Q)M 2 < N2~AP ~ 120' 



by the definition of N 2 . 



3.13 Switchings of colour 15 

Switchings of colour 15 are used to reduce the number of non-loop edges with multiplicity 3 
until this number is at most \N^/2\ . They are applied to multigraphs Q such that switchings 
of colours 1, . . . , 14 do not apply and e%{Q) > 

Let Q be such a multigraph. Then a switching of colour 15 is defined by a sequence of 
distinct vertices 

(t> , w , v 1 ,w 1 , v 2 , w 2 , t> 3 , w 3 ) 

such that 

• {^0,^0} has multiplicity 3 in Q, 

• each of {vj, Wj} is a simple edge in Q, for j = 1, 2, 3, 

• the edges {vq, Vj} and {wq, Wj} are absent in Q for j = 1, 2, 3. 
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The switching removes the edges {vj, Wj} for j = 0, 1, 2, 3 (setting the multiplicity of each to 
zero) and inserts the simple edges {i>o, Vj}, {wo, Wj} for j = 1,2,3, as illustrated in Figure [8 
for the case that D = 3. 

Suppose that the multigraph R can be produced by performing a switching of colour 15 
from Q. The number of ways to perform a switching of colour 15 in Q is at least 

2e 3 (Q) ([2e x (Q)] 3 - 0(^ ax ei(Q) 2 )) > e 3 (Q)M 3 , 



using Lemma 3T while the number of ways that R could be produced using a switching 
of colour 15 is at most Mf. Therefore we can take ai 5 (Q) = e 3 (Q)M 3 and b 15 (R) = Mf, 
leading to 

Mf 2Mf 1 

tt(g ' } ~ e 3 (Q)M 3 " iVaii^ " 120' 

by the definition of N3. 



3.14 Analysis of the switchings 



We now explain how to apply Theorem 2.1 to analyse these switchings on Q . (See the 



statement of Theorem 2.1 for the necessary notation.) 

We now define a directed graph r with vertex set V(r) = {vq : Q G Go}, where Vq is 
associated with the set S(vq) = {Q} containing one object. These sets are certainly disjoint. 
By a slight abuse of notation we will identify Vq and Q from now on, and write Q for both. 
The edge set of r is defined as follows: there is an edge QR in r with colour c G {1, . . . , 15} 
if and only if R can be obtained from Q using a switching of colour c. Since a switching is 
used only where no switching with a lower-labelled colour applies, there is at most one edge 
from Q to R in r. Hence the endvertices of an edge uniquely determine the colour of the 
edge. We take \ c {v) = je for all c, v, and so for each edge QR in r we have 

a{QR)= a{QR) = 15 a(QR). 
Mqr){R) 



Let Y, Z be as defined in (3.2), (3.3) and let X = (Jo- 



Lemma 3.2. With notation as established above, we have 

G(k,J,J*) = {l + O(kiJM))\g -Y\. 
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Proof. It follows from the analysis of the previous sections that a(QR) < and therefore 
a(QR) < for all edges QR in r. Moreover, we have provided at least one switching which 
can be performed from Q, for each graph Q in Qq — Z. Hence X, Y, Z satisfy the requirements 
of Theorem 12.11 and we can conclude that 



Y < 



i-d(y,y)' 



where 



a(W,W')= max a(P) 
Pev(w,w) 



for all W,W C Q . 

Since V(Y, Y) by definition has only non-trivial paths, and a(QR) < | for all edges in r, 
we know that a(Y,Y) < |, and therefore 

\Y\<2a(Y,Z)\G \. 

Now let P be a path in r from some Q G F to some element of Z, such that all internal 
vertices of P belong to Qq — (Y U Z). Let Qi? be the first edge in P. Then QR cannot have 
colour in {1, ... , 8}, since these switchings only produce graphs in Y. If QR is coloured with 
a colour in {9, 10, 11, 12} then a(QR) = 0(k 3 mSLX /M), so a(P) = 0(^ ax /M). 

The remaining possibility is that QR has colour 12 +r for some r e {1,2,3}, in which case 
P must contain at least [N r /2\ edges of colour 12 + r. Our analysis showed that a(e) < 
and hence that d(e) < |, for all such edges e. By definition [N r /2\ > |(logM — 1) for 
r G {1,2,3}. Therefore, 

a(P) < g-(^M-i)/2 = 0(M~ 1 ), 
which implies that a(Y,Z) = 0(/c^ ax /M). We conclude that 

\g -Y\ = (i + o(k 3 m jM))\g \. 

But g - Y C ^(Jfe, J, J*) C ^o, and hence 

G(k,j,r) = {i + o(kijM))\g -Yi 

completing the proof. □ 
It remains to obtain an asymptotic expression for g$ — Y, which we do in the next section. 
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4 From pairings to multigraphs 



In this section we work in the pairing model (also called configuration model), which we 
now describe. This model is standard for working with random graphs of fixed degrees: see 
for example [7]. Consider a set of M points partitioned into cells Ci, . . . , c n , where cell q 
contains hi points for % — 1, 2, . . . , n. Take a partition P (called a pairing) of the M points 
into M/2 pairs with each pair having the form {y, z} where y G c, and z G Cj for some 
The set of all such pairings, of which there are M\/(M/2)\ 2 M//2 , will be denoted by C(k). 

A loop is a pair whose two points lie in the same cell, while a link is a pair involving 
two distinct cells. Two pairs are parallel if they involve the same cells. A parallel class is a 
maximal set of mutually parallel pairs. The multiplicity of a parallel class (and of the pairs 
in the class) is the cardinality of the class. As important special simple pair is a 

parallel class of multiplicity one, a double pair is a parallel class of multiplicity two and a 
triple pair is a parallel class of multiplicity three. 

Each pairing gives rise to a multigraph in Q{k) by replacing each cell by a vertex, and 
letting the multiplicity of the edge {v , w} equal the multiplicity of the parallel class between 
the corresponding cells. 

Let Ci t d,t be the set of all pairings in C(k) with exactly I simple loops, exactly d double 
pairs and exactly t triple pairs, but with no loops of multiplicity greater than one and no 
links of multiplicity greater than three. If G G Go — Y can be formed from a pairing P G Cej,t, 
then exactly 

n 

pairings in C(k) give rise to G. Now defining 

w(£,d,t) = 2 i+d 6 t \C iAt \, 

we can write 

^ i=l ' i=0 d=0 t=0 

Hence it suffices to obtain an asymptotic expression for the above sum. 
We will need the following two summation lemmas proved in [5]: 

Lemma 4.1 ([21 Corollary 4.3]). Let < Ai < A 2 and P>i < B 2 be real numbers. Suppose 
that there exist integers N, K with N > 2 and < K < N, and a real number c > 2e such 
that < Ac < N - K + 1 and \BN\ < 1 for all A e [A U A 2 ] and B G \B U B 2 ). Further 
suppose that there are real numbers 5i, for 1 < i < N , and ji > 0, for < i < K, such that 
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E}=iN<E^7j[«1i< J/orl<<<JV. 

Given A(i), . . . , A(N) G [Ai,A 2 ] and B(l), . . . , B(N) G [B 1 ,B 2 ], define n ,ni,...,n N by 
uq = 1 and 

^ = ^(l-( J -l)5( ! ))(l + i i ) 

ni-i * 

/or 1 < « < iV ; mf/i the following interpretation: if A(i) = then nj = for i < j < N . 
Then 

N 



Ei<^2,ni< E 2 , 



where 



i=0 



K 

27 a = exp^x - 1^2-4^7,(3^)^) - \{2e/c) N 

j=0 

K 

U 2 = exp(A 2 - \A\B X + \A\B\ + 4^ 7i (3A 2 )') + |(2e/c)*. □ 



j'=o 

Lemma 4.2 ([51 Corollary 4.5]). Let N > 2 be an integer and, for 1 < i < N, let real 
numbers A(i), B(i) be given such that A(i) > and 1 — (z — l)-B(i) > 0. Define A\ = 
minf =1 A(i), A 2 = max^ A{i), C\ = mmf =1 A{i)B{i) and C 2 = m.ax.f =1 A{i)B{i). Suppose 
that there exists a real number c with < c < | such that max{v4/iV, \C\} < c for all 
A G [A\, A 2 ], C G [Ci, C 2 ]. Define n , . . . , by no — 1 and 

for 1 < i < N, with the following interpretation: if A(i) = or 1 — (i — l)B(i) = 0, then 
nj = for i < j < N . Then 

N 
i=0 

where 

E x = exp(A 1 - \A X C 2 ) - (2ec)* 

E 2 = exp(A 2 - \A 2 C X + \A 2 Cl) + (2ec) N . □ 

First we perform a summation over the number of edges of multiplicity three. 
Lemma 4.3. Uniformly for < d < N 2 and < I < N\, we have 

N z ...re a +\ / m2 



w(i,d,t) _ J Mj \ 
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Proof. We will apply Lemma 4.2 Let t' be the first value of t < N 3 for which Cg d,t — 0, or 



t' — N 3 + 1 if there is no such value. By [T2l Lemma 4.1] we have 

w{£, d, t) _ 6 \Ct At \ _ Ml - a t k 2 max {k 2 max + £ + d + t- l)M 3 



w(£,d,t-l) \C tAt -i\ 2tM 3 



(4.2) 



where ot t = a t (£,d) is bounded independently of £, d and M, for 1 < t < t' . Note that our 
expression for N 3 is no smaller than, and at most a constant factor larger than, that used 
in [12]. Examination of the proof of [121 Lemma 4.1] shows that the expression given there 
for \Ci )d>t \/\Ct,d,t-i\ also holds for 1 < t < N 3 . 
For 1 < t < t', define 

A(A _ Ml-a t k 2 milx (k 2 max + £ + d)M 3 _ a t k 2 max M 3 



2M 3 ' v ' 2M 3 

Now M 3 > since 1 < t < f. If a t < then by definition A(t) > M|/2M 3 > 0. Now 



suppose that a t > 0. Then C{t) > 0. Note that the left hand side of (4.2) is positive, and 
multiplying by t we see that A(t) — (t — l)C(t) is positive, which implies that A(t) > 0. 
Therefore A(t) > whenever 1 < t < f. Define B(t) = C(t)/A(t) for 1 < t < t'. Also define 
A(t) = B(i) = for if < t < N 3 . 

Define Ai,A2,Ci,C<i by taking the minimum and maximum of A(t) and C(t) over 1 < 

Let A G L4i,A 2 ] and C G [Ci, C2], and set c = ^. Since 



t < N 3 , as in Lemma 



4.2 



A = MI/2M 3 + o(l) and C = o(l), we have that m&x{A/N 3 , \C\} < c for M sufficiently 
large, by the definition of N 3 . 

Therefore Lemma 4.2| applies and says that 



§ = BXP & + °«»(^ + < + ")/ M2 )) + O((e/40)*)- 

Since £ + d < N\ + N 2 = 0(/c^ ax + log M), the relative error term is 0(k^ aax /M), as claimed. 
Finally, (e/40) 3 < (e/40) logM < M~ 2 . Since the sum we are estimating is at least equal 
to one, the additive error term is covered by the error term inside the exponential. This 
completes the proof. □ 

Next we perform a summation over the number of simple loops. 

Lemma 4.4. Uniformly for < d < N 2 , we have 

A w(£,d,0) _(M 2 fk max d h*\\ 
^^(0,d,0)- eXP U V M + M J J 
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Proof. Let £' be the first value of £ < N\ for which Cg^o = 0, or = A^ + 1 if there is no 
such value. By [T2], Lemma 4.2], for 1 < £ < £' we have 

w(£, d, 0) 2 \C W \ M 2 - (3 e (k 3 m , x + k max d + k 2 max (£ - 1)) 



w(£-l,d,0) |C,_ Mi0 | £M 



(4.3) 



where = 0e(d) is uniformly bounded (independently of d and M). As in Lemma 4.3 above, 
our expression for N± is no smaller than, and at most a constant factor larger than, that 
used in [12J. Examination of the proof of [121 Lemma 4.2] shows that the expression given 
there for |C^<z,o|/|C^-i,<i,o| also holds for 1 < £ < Ni. The proof is completed by adapting the 
argument used for Lemma 4.3, above. 

For 1 < £ < £' , define 

U M ' M ' 

By the same argument as in the previous lemma, A(£) > for £ < £'. Define B(£) = 
C(£)/A(£) for 1 < £ < £' . Also define A(£) = B(£) = for £' < £ < N u and c = i. 
Lemma 4.2| now applies to give the desired result. □ 

Next we perform a summation over the number of edges of multiplicity two. The expo- 



nential factor in the summand corresponds to the error term from Lemma [4T4] which depends 
on d. 

Lemma 4.5. For any constant p have 

A w(0, d, 0) fpk max d\ _ ( Ml Mj \ 
|^(0,0,0) eXI \ M )- eXP {2M 2 2M* +U[k ™ JM) )- 

Proof. We divide the proof into two cases, following the division used in [T2"j Lemmas 4.3- 
4.4]. Again, our expression for N 2 is no smaller than, and at most a constant factor larger 
than, that used in [12] . Examination of the proofs of [T2l Lemmas 4.3-4.4] shows that the 
expression given there for |Co,d,o|/|Co,d-i,o| also holds for 1 < d < N 2 . 



For < d < N 2 , define 



w(0,d,0) f pk max d 



m d = 77, n eX P 



w(0,0,0) M 



First suppose that M 2 < M. From [T2"j Lemma 4.3] we know that for 1 < d < N 2 with 
C , dfi ^ 0, 

m d _ 2 exp(pk max /M) \C 0A0 \ _ M| ( A; max (^max + d ) M 2 



md-i \Co,d-i,o\ 2dM 2 \ dM* 



25 



The current lemma follows from this, arguing as in the proofs of the previous two lem- 
mas. (Note that the term M|/2M 3 in the answer is absorbed into the error term, by the 
assumption that M 2 < M.) 

Now suppose that M 2 > M. By [TJl Lemma 4.4] we have for 1 < d < N 2 with C ,d,o 7^ 0, 
m d 2 exp(p/c max /M) \C 0A0 \ A(d) 



m<i-i \Co,d-ifi\ d 



(l-(d-l)B) (1 + 5, 



d 



where 



Ml 2M^M 3 Ml M 2 4 

A ^ = 2M! + T - 2i^ - W + 

8 16M 3 

M + M| ' 

5 d = 0((d - 1) /M 2 ) uniformly for 1 < d < N 2 . 

(Note that A(d) exp(pk nmx /M) equals A(d) within the precision afforded by the error term, 
since k max A(d)/M = 0(/c max M|/M 3 ) = 0(fc max /M).) Again, our N 2 is no smaller than, 
and at most a constant factor larger than, the value of N 2 used in [12], and examination 
of the proof of [12], Lemma 4.4] shows that the expression for |Co,d,o|/|Co,d-i,o| also holds for 
1 < d < N 2 . 

Let d' be the first value of d < N 2 for which Co,d,o = 0, or <i' = iV 2 + 1 if no such value 
of d exists. Define n d = m^jm^-x for < d < d' and n d = for d > d'. We will apply 



Lemma 4.1 to ^2 d nd. 

Let A\ = min^ A(d) and A 2 = maxdA(d), with the minimum and maximum taken over 
1 < d < N 2 , and let B\ = B 2 = B. Define c = 80 and K = 2. Then the conditions of 



Lemma 4.1 apply with 70 = 7i = and some value of 72 which satisfies 72 = 0(1/M 2 ). 
Application of that lemma now gives the desired result. □ 

Recall the definition of yi, x 2 , £3 G {0, 1} given just before the statement of Theorem |1.1[ 
By combining the last three summations we obtain the following. 

Lemma 4.6. As M — > 00 and fc max = o(M), we have 

aft ^} ^ d f) , M M 2 M 2 \ 

d=0 £=0 t=0 v ' ' y v 7 

Proof. From Lemmas 43] and 4.4, we have 



^ ^ w(£,d,t) _ w(0,d,0) / M 2 Mj fk^d k^ 
w(0,0,0) w(0, 0, 0) \ M 2M 3 \ M M 



=0 t=o 
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uniformly over d. 

If X2 = 0, we are finished. If x 2 = 1, choose two constants pi, pi such that the actual value 



of the error term 0(k mSLX d/M) in (4.4) lies in [pik max d / M , P2k max d/M] for all d G {1, . . . , N 2 }. 
Applying Lemma [4^5 with p G {pi, P2] and noting that the result does not depend on p within 
the precision given by the error term, we are done. □ 

Finally we can prove our main result. 



Proof of Theorem \l.l\ From Lemma 3.2 and Lemma 4.6 together with (4.1), we obtain 



G(k, J, J*) = (1 + O(kljM)) |0o - Y\ 



ki\ ■ ■ ■ k n \ 

1 n d=o e=o t=o 



w(0,0,0) / M 2 Mf . . Mf n/t3 



Now w(0, 0, 0) = |Co,o,o|; an d h follows from [12], Lemma 5.1] that 

Ml f_M 2 _ Ml _M 2 2 M 3 Ml Ml 3 

1 °' ' ' ~ (M/2)!2^/2 GXP \2M ~UP~ ~2M^~ + 4AP + GAP + °^ JM) 

Combining these two expressions completes the proof. □ 



5 Comparison to a naive model 

Let p G (0, 1) be a probability which we will define later. Abusing notation slightly, we 
define the functions 

jeJ jeJ* 
and the probability generating functions 

m .M, 9(z) 



J(l)' av ' J*(l)' 

Consider a random symmetric nonnegative integer matrix A = (a^-) which is created as 
follows: 

• All entries on and above the diagonal are independent; 

• For i — 1, . . . , n, let an be a randomly chosen element of J*, with Pr(a„ = b) = p b /J*(l) 
for all b G J*; 
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• For 1 < i < j < n let randomly chosen element of J, with Pr(ajj = b) = p b /J(l) 
for all b G J; 

• Finally, let = for 1 < % < j < n. 

We write Pr(-) to mean probabilities generated by the above procedure. If Aq is a fixed matrix 
which corresponds to an element of Q(k, J, J*) then the probability that A is produced by 
the above procedure is 

M/2 

Pr(^o) = ^ ■ (5.1) 

This follows as the probability of Aq depends only on the sum of the diagonal and above- 
diagonal entries (that is, on the number of edges of the corresponding multigraph). 

For i = 1, . . . , n, let £i denote the event that the sum of row i equals ki, recalling that 
diagonal entries are weighted by a factor of 2 as in ( 1.1 ). Now 

r( u j m _ Pr(£iA£ 2 A---Ag n ) 

since Pr(-) is uniform on (the set of matrices corresponding to) Q(k, J, J*). By (incorrectly) 
assuming that the events £\, . . . ,£ n are independent, we obtain the following naive estimate 
of G(k, J, J*), using Q : 

niLiPr(^) 



G p (k, J, J*) 

J(l)(S)j*(i)» ™ r 2N 



y i=l 
n 

-Ml"}. -r/, S-( n \ ' 

V 



M ' 2 J(l)-U) Jl [z kl ] {J{z) n - l r{z 2 )). (5.2) 



i=l 



Let k = M/n denote the target average row sum. The value po of p which makes the 
expected row sum equal to k satisfies 



k (1 — 2x 2 )k 2 

Po=~ + ~ + 0(k n 2 ). 5.3 

n n z 



We write G nSu i vc (k, J, J*) for G Po (k, J, J*) from now on. (As we shall see, the exact value of 
the 0(k/n 2 ) term does not affect Theorem 5.1 below, within the stated error bound, even 
though some of the intermediate formulae we give in the proof are affected.) 

We now compare the expression of Theorem 1 1 . 1 1 with the naive estimate G nsiiive (k, J, J*). 
For future applications, it will be convenient to express the result in terms of the scaled 
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central moments /i 2 , fi 3 defined by 



* M 

i=l 



n 

,7 5> " *)' 



for r = 2, 3. 

Theorem 5.1. Under the conditions of Theorem \l.l\ we have 

G(k,J, J*) = V2G naivc (k, J, J*) exp f\(l - /i 2 ) (1 + 2x 2 + /i 2 (l - 2x 2 )) 



6/i 2 /i 3 (x 3 - x 2 ) - (4 + 3^2(^2 - 2 /i 3 ) + 2/13(3x3 - 3x 2 + 1) 



+ ^(9*3 - 9x 2 - 1) + Q{kljM) ) . (5.4) 



2n 12M 

2n 2 

Proof. Suppose that 

k (1 - 2x 2 )k 2 ck 
P=~ + - 2^— + ^ 

where c = c(n) = 0(1). By direct computation we find that 

P M/2 = ( ~) M/2 ew( \k\l - 2x 2 ) ~^ + C 4 + O(kiJM) ) (5.5) 

and 



nj 4n ' 2 



J(l)l» = exp( r)p - (1 ~ 2x ^ p2n2 + (1 - 3X2 + 3X3y " 2 + 0(^ ax /M) 



2/4 6 
= exp(lM - \k{2x 2 k + 2 - k) + ^ - k ^+^M + o^M)) . (5.6) 

To calculate the product in ( |5.2 ) we consider two cases, depending on the value of x 2 - 
First assume that x 2 = (that is, the value 2 is not permitted off the diagonal). Then 

J(z) n ~ l r(z 2 ) = (l+pz) n ~ l H{z), 

where 



Consequently, for any integer r with < r < /c max we have 
[z r ] JizT- 1 J*(z 2 ) = ^(11 % r ~ j WW 



i=Q 

fn - 

r ) r \r — 2/ i7± ^ Vr — 3/ ° v ^ V r 



V + n )yiP + Wn-l)p r + P r A (5-7) 
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where 

r / -i\ / n \-l 



4=e(;: j 1 )( b ; 1 )V' m * w - 



2\-l 



i=4 

Now define 

Notice that the coefficients of H + (z), when expanded as a Taylor series in 2, are nonnegative 
and dominate those of H(z). Also notice that (™Zj) ("^) 1 = 0(l)(r/n)^ uniformly for 
4 < j < r, since r = c^n 1 ' 2 ). Therefore, setting a = r/pn, 



\A\ = 0(l)J2<* j [z j ]H + (z) 



J=4 

= O (# + (a) - 1 - pa 2 - (n - l)pW) 
= 0{rklJM% 

where the last line follows since H + (z) — 1 — pz 2 — (n — l)p 3 z 3 has three zero derivatives at 
z = and its fourth derivative is 0{p 2 ) for 2 = 0(A; max ). Applying this bound to (5.7), we 
have 

\z T ] j{z) n ~ l r{z 2 ) 



r J \ kn n 2 

= «J exp (-r± + V + MiH^il + + 0(r ^ /M2) , . (5 . g) 

Now we assume that X2 — 1 and repeat the calculation. We have 

(\ 71—1 
(i-^)5> P v (1 + 5] 
i>0 ' i>l 

Arguing as before, we find that for any r with < r < fc maX ; 

[z r ] J(z) n -V*(z 2 ) 

n + r-2\ fn + r-4\ x fn + r - 5\ , 

V+l r _ 2 J2/1P +( r _ 3 l(a;3-l)(n-l)p 



r 

n + r-2\ (y x r{r-l) (x 3 - l)r 3 3 „, 

p exp + g 1f ~ 

r I \ kn n 



+ [ n + [ 2 ) P r O(rkiJM 2 ) 



f I 
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Combining the cases (5.8) and (5.9) we have, for i = 1, . . . ,n 

l + 2x 2 , Hi 



z ki \ J(z) n - X r{z 2 ) 



h le- 
ft p 



exp 



2n 



kn 



1 - 2x 



2n 



2 Vi_ 
kn 



1 + 6 ^ 2 ~ 6X3 ^ 3 + Q(^/M 2 )). 



Multiplying these n equations together gives 



n M p M 



exp I 



(M 2 {2x 2 - 1) -2M yxMa 



2n kn 
(6x 2 - 6a; 3 + 1)M 3 
6n 2 



(5.10) 



We now substitute (5.5), (5.6) and (5.10) into (5.2). Writing the result in terms of fx 2 , 
U3, M and n, using the identities 

M 2 = Mu 2 + (k - 1)M, 

M 3 = Mu 3 + 3(k - 1)M a 2 + (k- l)(Jfe - 2)M 

we find that 

C^naive (k, J, J ) 

1 /M\ M/2 / (2 yi + 2x 2 (/i 2 -l)-/i 2 )fc (2x 2 -l)P 
exp -2/1(1-^2) + + - 



4 



A;(6x 2 - 6x3 + l)(2/i 3 + 6/i 2 /c + k 2 
12n 



+ 0(C/M) . 



(At this point we can observe that no remaining terms depend on c, which verifies our earlier 
claim that the exact value of the 0(k/n 2 ) term in ( |5.3| ) does not affect the statement of this 
theorem.) The proof is completed by applying Theorem |1.1| and using Stirling's formula. □ 



.2 

u max 



Corollary 5.2. Under the conditions of Theorem 5.1, if u 2 = 0(M l l & ) then 
G{k,J,J*) 

= V2G naiye {k,J,J*) exp(i(l-0 2 )(l + 2z 2 +/z 2 (l-2x 2 ))+0(fe 

1 /2 

If the stronger bound \i 2 = 0(fc ma x) holds then 
G(k,J, J*) 

= V2G^ ivc (k,J,J*) exp(i(l-// 2 )(l + 2x 2 + M2 (l-2x 2 )) +0(klJM) 
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Finally, if k = (k,k,k . . . ,k) is a regular degree sequence with kn even, then 



G(k, J, J*) = s/2 exp(|(l + 2x 2 ) + 0(k 2 /n)) G naivc (fc, J, J*). (5.11) 
Proof. For the first two statements, we just need to check that all additional terms inside 



the exponential factor in Theorem 5.1 are covered by the claimed error bounds. For this 



it is useful to note that /x 2 < k max and \ps\ < k max p, 2 . Finally, when k is a regular degree 



sequence we have /x 2 = ^3 = and (5.11) follows. □ 



The constant ±(1 + 2x 2 ) in (jSJJj) was previously noted for simple sparse regular graphs 
in [12] , and in [3] for simple sparse regular graphs with loops allowed. Interestingly, with 
different error terms, the same constant was observed in [12] for dense simple regular graphs, 
in [3] for dense simple regular graphs with loops allowed, and in [9] for dense symmetric 



integer matrices with zero diagonal. In these three cases, we conjectured that (5.11) holds 
(with some vanishing error term) for all degree sequences except for the two extreme cases 
of graphs with no edges and graphs with all possible edges. The corresponding conjectures 
may be less likely to be true in full generality for all J, J*. 



From the first expression of Corollary 5.2, we see that G{k, J, J*) is closely approximated 
by y/2 G na _[ ve (k, J, J*) whenever /i 2 is close to 1. It appears that, in the random matrix model 
described at the start of this section, /i 2 will be concentrated around 1 with high probability 
whenever p tends to zero slowly enough to ensure that M — > oo with high probability, but 
quickly enough to ensure that fc^ ax = o(M) with high probability. If so, this will lead to a 
model for the degree sequences of sparse multigraphs analogous to that obtained by McKay 
and Wormald [13] for graphs, and by McKay and Skerman [TU] for bipartite graphs and 
directed graphs. Details will be given in a future paper. 
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